Degenerate Fermi Gas of ^''Sr 
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We report quantum degeneracy in a gas of ultra-cold fermionic *^Sr atoms. By evaporatively 
cooling a mixture of spin states in an optical dipole trap for 10.5 s, we obtain samples well into 
the degenerate regime with T/Tf = 0.26l'^Qg. The main signature of degeneracy is a change in the 
momentum distribution as measured by time-of-flight imaging, and we also observe a decrease in 
evaporation efficiency below T/Tf ~ 0.5. 



Quantum degeneracy of ferniions in dilute atomic gases 
[l| is currently one of the most active areas of physics 
research. Tunable interactions [3| have allowed an explo- 
ration of the BEC-BCS crossover [3| as well as study of 
phenomena in the unitary regime such as superfluidity 
in systems with [J, |5[ and without [Q^ spin imbalance. 
Research on Fermi gases in optical lattices [9| can make 
direct connections to the properties of electrons in solids 
and realize important models like the Hubbard Hamilto- 
nian [10], but with new capabilities for controlling sys- 
tem parameters such as density, interaction strength, and 
dimensionality. Current searches are underway for ana- 
logues of high temperature superconductivity |llj and 
spontaneous Neel magnetism [12l |. 

Here we report quantum degeneracy of the fermionic 
isotope of strontium, ^''Sr. Five other species of fermions 
have been brought into the quantum degenerate regime 



C^^K [131, ^Li [IJ, 3He [15|, ^^^Yb and ^^^Yb [161), but 



^'^Sr has several properties which make quantum degener- 
ate samples of this type particularly interesting. *^Sr has 
a very large nuclear spin, 1=9/2, which may allow stud- 
ies of novel magnetic phenomena due to enlarged SU(N) 
symmetry of the interaction Hamiltonian for A^ = 2/ -I- 1 
17h19I|. High resolution spectroscopy technologies are 
the most advanced in strontium because of the use of 
narrow intercombination transitions in ®^Sr for optical 
frequency standards [201 > ^'^'i these tools have motivated 
proposals for applications in quantum information [21- 
l23l and quantum simulation of many-body phenomena 
[171 l24l . Strontium also has a number of stable bosons 
which have recently been brought into the quantum de- 
generate regime [25- 27], which makes Fermi-Bose mix- 
tures with relatively small mass differences available [28| . 
There is also the potential for manipulating interactions 
on small spatial and t empora l scales with low-loss optical 
Feshbach resonances 



29,3 



Details about our apparatus can be found in [26 . 



3ll . |32]- Atoms are trapped from a Zeeman slowed 
beam in a magneto-optical trap (MOT) operating on the 
(5s^)"'^iS'o — (5s5p)^Pi transition at 461 nm. Since this 
transition is not closed, approximately 1 in 10^ excita- 
tions results in an atom decaying through the ipsbd)^ D2 
state to the {hsbpf'P2 state, which has an 9 minute life- 
time [33! and can be trapped in the quadrupole mag- 
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FIG. 1. (color online) Partial level diagram for Sr (- -) and 
*^Sr ( — ) including hyperfine structure and isotope shifts. For 
*^Sr, total quantum number F is indicated. 



netic field of the MOT [3J-|37|. The magnetic trap has 



a lifetime of about 25 s, which is limited by background 
pressure and blackbody radiation [36|. This allows us to 
accumulate atoms over a period of 30 s and trap a sig- 
nificant number in spite of the low natural abundance of 
«"Sr (7%). 

After accumulation, atoms in the ^P2 state are 
returned to the ground state via excitation to the 
(5s4d)^£'2 state [32|. This repumping is achieved by ap- 
plying 3W/cm^ of 3.012 /im light for 60 ms. The tran- 
sition has hyperfine structure due to the nuclear spin of 
/ — 9/2, which spreads the transition over ~3GHz. In- 
dividual hyperfine transitions are not resolved in the re- 
pumping efficiency curve because of the high intensity of 
the repump laser and len gth of time over which the re- 
pumping laser is applied 38| . The 3 /im laser is tuned 
1 GHz blue of the ^^Sr centroid, which effectively re- 
pumps transitions from the F — 11/2 and F = 13/2 ^P2 
levels. During the repumping stage, the 461 nm MOT is 
left on so that atoms returned to the ground state are 
recaptured and cooled. We typically recapture 3 x 10^ 
atoms at temperatures of a few milliKelvin. It should 
be possible to improve this number by modulating the 
repumping laser frequency to excite atoms in all of the 
F levels. 



Following this step, the 461 nm light is extinguished 
and 689 nm light is applied to drive the {bs'^^So — 
{5s5p)^Pi transition and create an intercombination line 
MOT. Similar to [3^, two frequencies are applied (Fig. 
[T]). The two frequencies are required because the Zee- 
man shift is much larger in the excited state than the 
ground state, causing some ground-state mp levels to be 
anti-trapped. A "trap laser", slightly red-detuned from 
the F ~ 9/2 -^ F = 11/2 hyperfine transition, cre- 
ates a stable trapping force on average for atoms only 
in the presence of a "stir laser" red-detuned from the 
F = 9/2 -^ F — 9/2 transition. The stir laser effectively 
randomizes the populations because the Zeeman shift for 
this transition is much smaller than for the trapping tran- 
sition. For a strong transition, such as the ^Sq—^Pi tran- 
sition at 461 nm, this mixing is accomplished by a single 
trapping laser, but an extra laser is required for the in- 
tercombination transition because excitation is relatively 
slow. 

To increase the initial capture rate from the 461 nm 
MOT, both the stir laser and trap laser are artificially 
broadened so that they are resonant with a wider veloc- 
ity range of atoms. This is achieved by modulating the 
frequencies of the acousto-optic modulators generating 
the beams at a frequency of 30 kHz with a peak-to-peak 
amplitude of 660 kHz. Initial detunings are -1.2 MHz 
and -900 kHz from the stir and trap transitions respec- 
tively, and intensities are ImW/cm^ and 2.2mW/cm^. 
This allows trapping of up to 70% of the atoms in 
the intercombination-line MOT. During a 400 ms cool- 
ing time, the detuning, dither, and intensity are reduced 
to -90 kHz, 200kHz, and O.SmW/cm^ for the trap laser 
and -600 kHz, 200 kHz and O.gmW/cm^ for the stir laser. 
Simultaneously, the magnetic field gradient is increased 
from 0.2G/cm to 1.9G/cm. This increases spatial con- 
finement and transfer to the optical dipole trap (ODT). 

After the cooling and compression stage, an ODT, 
consisting of two crossed laser beams, is overlapped for 
100 ms with the intercombination-line MOT with a mod- 
est power (3.9 W) per beam. A single beam derived from 
a 21 W multimode, 1.06 fiia fiber laser is recycled through 
the chamber to form the ODT. The resulting trap has 
equipotentials that are nearly oblate spheroids, with the 
tight axis close to vertical. Within the trapping region, 
each beam has a waist of approximately 90/L(,m. Dur- 
ing this loading time, the dither amplitudes are reduced 
to zero, and the intensities and detunings of the MOT 
beams are further reduced to approximately 3^W/cm^ 
and -30 kHz. 

After the extinction of the 689 nm light, the ODT 
power is ramped in 30 ms to 7.5 W per beam to obtain 
a trap depth of 25 /iK. At this point, we typically trap 
3 X 10^ atoms at 7//K with a density of 2.5 x 10^^ cm~^. 
No attempt is made to spin polarize the sample. Ex- 
periments using ^^Sr for optical frequency standards [40| 
found that the intercombination-line MOT produced a 



sample with roughly equal populations of ground-state 
magnetic sublevels. For a conservative estimate of the 
phase-space density, we will assume equal populations in 
all levels. This corresponds to an average initial colli- 
sion rate of 200 s^^ and a phase-space density of approx- 
imately 10^-^. As discussed below, exact determination 
of the polarization is not necessary to establish Fermi 
degeneracy of our sample and will be the topic of future 
studies. The sample lifetime in a static ODT, presum- 
ably limited by background gas collisions, is 30 s, which 
is sufficient for evaporative cooling. 

After loading into the ODT, we begin forced evapora- 
tion by lowering the power in the ODT according to the 
formula P = Po/{l + t/T)^ + PoSsot with time denoted by 
L_B = 1.4, and r = 1.5 s. This trajectory was designed 
41| without Poffsct to yield efficient evaporation when the 
effect of gravity can be neglected. Gravity is a significant 
effect in this trap for Sr, and to avoid decreasing the trap 
depth too quickly at the end of the evaporation, we set 
-Poffsct — 0.7 W which corresponds to the power at which 
gravity causes the trap depth to be close to zero. 

For diagnostics, the ODT laser is extinguished, and 
atoms are allowed to expand for a time-of-fiight (TOF) 
of between 20 and 22 ms. We then obtain absorption im- 
ages using a linearly polarized laser on resonance with the 
^So{F = 9/2) - iPi(F = 11/2) transition. The intensity 
is 0.02/s, where Is = 42mW/cm^ is the saturation in- 
tensity, and the exposure time is 35 /is. Due to the small 
splitting between the F = 11/2 and F — 9/2 ^Pi states, 
the imaging beam excites transitions to both states. For 
a given density of atoms, the optical depth of the sam- 
ple depends upon the population distribution of ground 
state magnetic sublevels. Numerical solution of the op- 
tical pumping rate equations show that atoms quickly 
approach a steady state distribution. For the small ini- 
tial polarization in our samples the variation in optical 
depth is relatively small and an average absorption cross 
section of cr = 3A^/47r, where A = 461 nm, is accurate 
to within 10%. This diagnostic does not provide enough 
information to distinguish atoms in different magnetic 
sublevels. 

Figure [2] shows data and fits for one-dimensional slices 
through the center of TOF absorption images along axes 
perpendicular to the imaging beam and to each other. 
Position has been converted to velocity assuming ballistic 
expansion and long expansion time. Ballistic expansion 
is a good approximation because the mean field interac- 
tion energy is negligible compared to the kinetic energy, 
and the expansion times used are much longer than the 
relevant time scale, l/w, where / = a;/27r « 100 Hz is the 
typical mean trap oscillation frequency. 

Well outside the quantum degenerate regime, the 
Maxwell-Boltzmann and Fermi-Dirac momentum distri- 
butions are not distinguishable. To fit the full TOF im- 
ages accurately when the atoms are in the quantum de- 
generate regime, we must use the Fermi-Dirac expression 
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FIG. 2. (color online) Velocity distributions along axes per- 
pendicular to imaging beam. At i = 4.1s of forced evapora- 
tion (top), the data is fit well with a Maxwell-Boltzmann(- 
-) or Fermi-Dirac distribution (-). At 10.5 s, classical-particle 
statistics overestimates the population at low velocities, while 
Fermi-Dirac statistics for T/Tp = 0.26 accurately fit the data. 



for the two-dimensional momentum distribution [42 



T^{Px,Py) 






-Li2[-C,e' 



(1) 



where Li„ is the poly-logarithmic function of order n, 
Pxiy are momenta transverse to the imaging beam, and 
C, is the fugacity. The fugacity is related to the sam- 
ple temperature (T) and Fermi temperature {Tp) by 
Lij,[—Q] = ~ 6(T/T )^ ' ^^ ^^^^ method yields a measure 
of the quantity T /Tp without an explicit assumption on 
the spin polarization. Formally, the signal is a sum of 
distributions for each magnetic sublevel, where each may 
have different fugacities. Because of the weak depen- 
dence of Tp = {6Ny/^huj/kB on the number of atoms 
in a single quantum state (N), atoms in magnetic sub- 
levels with the smallest T /Tp values will dominate the 
signal. Sublevels with significantly lower population and 
higher values of T/Tp will contribute less to the signal, 
but they do cause the extracted value of T/Tp to be an 
upper bound of the value for the most populated state. 
The weak dependence of the velocity distribution on the 
fugacity far from degeneracy makes extracting T/Tp dif- 
ficult for hotter clouds. 

When fitting the TOF images using Maxwell- 
Boltzmann statistics, the entire image is fit with a 
Gaussian. This fit determines the number of atoms 
[42| and the size of the cloud, s, where the density is 
n oc cxp(— r^/2s^). To determine the temperature, the 
high velocity wings of the distribution, consisting only 
of atoms at radii greater than s, are fit to a Maxwell- 
Boltzmann velocity distribution. The deviation of the 
full distribution from a Gaussian is never large and the 



FIG. 3. (color online) Evaporation of Sr. T/Tf is obtained 
from either a fit of the velocity distribution to a single Fermi- 
Dirac distribution (DFG), or a determination of number and 
temperature from fits to a Maxwell-Boltzmann distribution 
and knowledge of the trap oscillation frequencies (MB). (In- 
set) Variation oiT/Tp with atom number. T/Tf is calculated 
using the latter of the two methods. Evaporation is efficient 
early in the trajectory and T/Tf decreases significantly as 
atoms are lost. But T/Tf plateaus near the limit of about 
~ 0.25 as the evaporation efficiency decreases. 



high velocity wings are less sensitive to the changes in the 
shape of the distribution as the atoms become degener- 
ate 13| . Numerical simulations of velocity distributions 
show that this method of thermometry is accurate to 
within 10% for T/Tp > 0.4 and becomes less accurate 



for more degenerate samples. The temperature is always 
overestimated. Combining this measure of temperature 
with knowledge of the ODT oscillation frequencies pro- 
vides a measurement of T/Tp. Here we assume equal 
population of the magnetic substates, which makes this 
an upper bound on T/Tp for the magnetic sublevel with 
the largest population. 

After 4.1s of forced evaporation, the sample has 
reached T/Tp = 1 as measured by both methods of fit- 
ting, although the uncertainty is large for the Fermi- 
Dirac fitting. Both the Maxwell-Boltzmann and the 
Fermi-Dirac fitting describe the data well (Fig. [2]). After 
10.5 s the Maxwell-Boltzman distribution overestimates 
the number of atoms with small velocity in the central re- 
gion of the cloud. However, the Fermi-Dirac fitting yields 
an excellent description of the data with a fugacity cor- 
responding to T/Tp = 0.26t'o6- The non-Gaussian char- 
acter of the distribution is a clear signature of the onset 
of quantum degeneracy and the limiting of occupancy of 
lower energy levels due to the Pauli exclusion principle. 

Figure [3] shows the evolution of T/Tp during the forced 
evaporation, determined by the two methods described 
above. The temperature starts at 4 /iK and reaches 30 nK 



after 10 s of evaporation. We achieve a limiting value 
T/Tp = 0.25 with this evaporation trajectory, and the 
agreement of the two methods of determining this quan- 
tity indicate that the assumption of equal ground state 
sublevel populations is a reasonable approximation for 
this data. 

Figure [3] (inset) provides information on why lower val- 
ues of T jTp are not observed. Evaporation initially pro- 
ceeds efficiently through collisions of ^'^Sr atoms in differ- 
ent magnetic sublevels. T/Tp drops by a factor of almost 
10 for a loss of about a factor of 10 in atom number. At 
values below T/Tp = 0.5, the evaporation becomes very 
inefficient. This signature has been seen in evaporation 
of 4°K [H] and "^Yb JJ] and ascribed to Pauh block- 
ing of collisions when atoms have reached the quantum 
degenerate regime. 

We have described the creation of a quantum degen- 
erate Fermi gas of *^Sr. This result opens possibilities 
of future studies involving quantum degenerate Bose- 
Fermi mixtures [26|, |27| , quantum degenerate Fermi gases 
with large ground state degeneracy 17H19J , and quantum 
computing with alkaline-earth metal atoms [2lM23|. 

Acknowledgements This research was supported by 
the Welch Foundation (C-1579), National Science Foun- 
dation (PHY-0855642), and the Keck Foundation. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 

[8] 

[9] 

[10] 

[11] 
[12] 

[13] 
[14] 

[15] 



S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. Mod. 

Phys. 80, 1215 (2008). 

C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 

Mod. Phys. 82, 1225 (2010). 

I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys.| 

80, 885 (2008). 

G. B. Partridge, W. Li, R. I. Kamar, Y.-a. Liao, and 

R. G. Hulet, Science 311, 503 (2006). 

M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and 

W. Ketterle, Science 311, 492 (2006). 

K. M. O'Hara, S. L. Hemmer, M. E. Gehm, S. R. 

Granada, and J. E. Thomas, Science 298, 2179 (2002). 

C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, 

S. Jochim, J. H. Denschlag, and R. Grimm, Science 305, 

1128 (2004). 

M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, C. H. 

Schunck, and W. Ketterle, Nature 435, 1047 (2005). 

M. Kohl, H. Moritz, T. Stoferle, K. Giinter, and 

T. Esslinger, Phys. Rev. Lett. 94, 080403 (2005). 

R. Jordens, N. Strohmaier, K. Gunter, H. Moritz, and 

T. Esslinger, Nature 455, 204 (2008). 

W. Hofstetter, J. L Cirac, P. ZoUer, E. Dernier, and M. D. 

Lukin, Phys. Rev. Lett. 89, 220407 (2002). 

L.-M. Duan, E. Demler, and M. D. Lukin, 

Phys. Rev. Lett. 91, 090402 (2003). 

B. DeMarco and D. S. Jin, Science 285, 1703 (1999). 

A. G. Truscott, K. E. Strecker, W. I. McAlexander, G. B. 

Partridge, and R. G. Hulet, Science 291, 2570 (2001). 

J. M. McNamara, T. Jeltes, A. S. Tychkov, W. Hoger- 

vorst, and W. Vassen, |Phys. Rev. Lett.| 97, 080404 

(2006). 



[16 

[17 

[18 
[19 
[20 
[21 
[22; 

[23; 
[24; 
[25; 

[26; 

[27 

[28; 

[29 
[30 
[31 

[32; 

[33 

[34; 
[35; 

[36 

[37 

[38; 
[39 
[40 

[41 
[42; 



T. Fukuhara, Y. Takasu, M. Kumakura, and Y. Taka- 

hashi, Phys. Rev. Lett. 98, 030401 (2007). 

A. V. Gorshkov, M. Hermele, V. Gurarie, C. Xu, P. S. 

Julienne, J. Ye, P. ZoUer, E. Demler, M. D. Lukin, and 

A. M. Rey, Nature Physics 6, 289 (2010). 

M. Hermele, V. Gurarie, and A. M. Rey, Phys. Rev. Lett. 

103, 135301 (2009). 

M. A. CazaUUa, A. F. Ho, and M. Ueda, New Journal of 

Physics 11, 103033 (2009). 

J. Ye, H. J. Kimble, and H. Katori, Science 320, 1734 

(2008). 

A. J. Daley, M. M. Boyd, J. Ye, and P. ZoUer, Phys. Rev. 

Lett. 101, 170504 (2008). 

A. V. Gorshkov, A. M. Rey, A. J. Daley, M. M. Boyd, 
J. Ye, P. ZoUer, and M. D. Lukin, Phys. Rev. Lett. 102, 
110503 (2009). 

I. Reichenbach, P. S. Julienne, and L H. Deutsch, Phys. 

Rev. A 80, 020701 (R) (2009). 

M. Foss-Feig, M. Hermele, and A. M. Rey, 

arXiv:0912.4762(2009). 

S. Stellmer, M. K. Tey, B. Huang, R. Grimm, and 

F. Schreck, Phys. Rev. Lett. 103, 200401 (2009). 

Y. N. Martinez de Escobar, P. G. Mickelson, M. Yan, 

B. J. DeSalvo, S. B. Nagel, and T. C. KiUian, Phys. Rev. 
Lett. 103, 200402 (2009). 

P. G. Mickelson, Y. N. Martinez de Escobar, M. Yan, 
B. J. DeSalvo, and T. C. KiUian, arXiv: 1003.3867(2010). 

A. Imambekov and E. Demler, |Phys. Rev. A| 73, 
021602(R) (2006). 

R. Ciurylo, E. Tiesinga, and P. S. Julienne, Phys. Rev. 

A 71, 030701(R) (2005). 

K. Enomoto, K. Kasa, M. Kitagawa, and Y. Takahashi, 

Phys. Rev. Lett. 101, 203201 (2008). 

Y. N. Martinez de Escobar, P. G. Mickelson, P. Pellegrini, 

S. B. Nagel, A. Traverso, M. Yan, R. Cote, and T. C. 

KiUian, Phys. Rev. A 78, 062708 (2008). 

P. G. Mickelson, Y. N. Martinez de Escobar, P. Anzel, 

B. J. DeSalvo, S. B. Nagel, A. J. Traverso, M. Yan, and 
T. C. KiUian, J. Phys. B 42, 235001 (2009). 

M. Yasuda and H. Katori, fhys. Rev. Lett., 92, 153004 
(2004). 

G. Ferrari, R. E. DrulUnger, N. Poll, F. Sorrentino, and 
G. M. Tino, Phys. Rev. A 73, 023408 (2006). 

S. B. Nagel, C. E. Simien, S. Laha, P. Gupta, V. S. 
Ashoka, and T. C. KiUian, Phys. Rev. A 67, 011401(R) 
(2003). 

X. Xu, T. H. Loftus, J. L. HaU, A. GaUagher, and J. Ye, 
J. Opt. Soc. B 20, 968 (2003). 

N. Poll, R. E. DruUinger, G. Ferrari, J. Leonard, F. Sor- 
rentino, and G. M. Tino, Phys. Rev. A 71, 061403(R) 
(2005). 

P. G. Mickelson, Trapping and Evaporation of ^"^ Sr and 
^^Sr Mixtures, Ph.D. thesis. Rice University (2010). 
T. Mukaiyama, H. Katori, T. Ido, Y. Li, and M. Kuwata- 
Gonokami, Phys. Rev. Lett. 90, 113002 (2003). 

A. D. Ludlow, The Strontium Optical Lattice Clock: Op- 
tical Spectroscopy with Sub-Hertz Accuracy, Ph.D. thesis, 
University of Colorado at Boulder (2008). 

K. M. O'Hara, M. E. Gehm, S. R. Granade, and J. E. 
Thomas, Phys. Rev. A 64, 051403(R) (2001). 

B. DeMarco, Quantum Behavior of an Atomic Termi 
Gas, Ph.D. thesis. University of Colorado (2001). 



